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ABSTRACT 

Motivated by recent observations of deuterated molecules in the interstel- 
lar medium, we examine the production of HD and D2 molecules on dust 
grain surfaces. A mechanism for the enhanced production of these deuterated 
molecules is studied. This mechanism applies on grain surfaces on which D 
atoms stick more strongly than H atoms, under conditions of low flux and 
within a suitable range of temperatures. It is shown that under these con- 
ditions the production rates of HD and D2 are greatly enhanced (vs. the H2 
production rate) compared with the expected rates based on the adsorption of 
gas-phase atomic abundances of D and H. The enhancement in the formation 
rate of HD is comparable with the enhancement due to gas-phase ion-molecule 
reactions in diffuse clouds. 
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1 INTRODUCTION 



The formation of molecular hydrogen on the surfaces of dust grains in the interstellar medium 
is a process of fundamental importance because the formation r ate in the gas phase can- 
not account for the observed abundance (|Gould fc Salpeteiill963h . Once produced on grain 
surfaces and ejected or desorbed into the gas, hydrogen molecules are a necessary compo- 
nent in the initiation of gas-phase chemical reaction networks that give rise to the chem- 
ical complexity observed in inte rstellar clouds. The formation of r nolecular hydrogen has 



been studied both theoreticallv (IHoUenbach &: Salpeter 



WimamTl968 : 



Smoluchowski 



] 198l[ 198.4 



1970 



Aronowitz fc Chang 



197 



1985 



HoUenbach et al 



Dulev fc WiUiam! 



1971 



1986 
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Pirronello & Avera 


1988; 


Sandford & AUamandola 


1993; 


Takahashi e 


; a 


1991 


3) and experi- 


mentallv ( 


Brackmann 


im 


Schutte et alJl97fil 


Ipirronello et al.lll997a 


hi 


999; 


Manico et al. 



I2nnif ). Although formation of H2 is clearly the most studied surface reaction pertaining to in- 
terstellar clouds, there are many other reactions that are thought to take place, for the most 
part involving H atoms and other species. Reactions are typically considered to occur via the 
Langmuir-Hinshelwood mechanism, in which weakly-bound surface species diffuse towards 
one an other. Some authors have considered different mechanisms, where chemisorption plays 



a role (^Cazaux fc Tielens 



The computational modeling of diffusive chemical reactions on dust grains in the in 



terstellar medium is typically done using rate equation models ( 


Pickles & Wi 


liams 


1977 




d'Hendecourt et al.ll985 


BrowJ 


1990 


Brown & Chan 


aleJl99n: 


Haseeawa et al. 


1992; 


Haseeawa & Herbst 


1993a 


y 


Casein et al.ll993l 


Willacv & WiUiamJl993l 


Shalabiea & Greenber Jll994 


) . As long 



as a sizable number of reactive species exists on a grain surface, the rate equation results are 
valid. However, in the limit of small average surface populations 1) of reactive species, 
rate equations are not always suitable, because they take into account only average den- 
sities and ignore the fiuctuatio ns as well as the discrete nature of the populations of the 



atomic and molecular species (iTielen 



Caselli et al, 



Shalabiea et al 



1995 



Tielens &: Hagen 



Stantcheva et al 



2m 



982; 



Charnlev et al 



1997 



20021 ). Depending on the 



parameters chosen, this limit can pertain mainly to small grains or can even to "classical" 
grains, namely grains of 0.1/i radius. For example, as the number of H atoms, A^h, on a 
grain fluctuates in the range of 0, 1 or 2, the H2 formation rate cannot be obtained from the 
average number alone. This can be easily understood, since recombination requires at least 
two H atoms simultaneously on the surface. 

Recently, a master equation approach was proposed that is suitable f or the simulation of 



diffusive chemica. 



Green et al 



\2m 



reactions on interstellar grains in the accretion limit (iBiham 



il .1120011 



Biham fc LipshtatI 120021 ). The master equation takes into account both 



the discrete nature of the reactive species as well as the fiuctuations in their populations. 



Green et al. 



( 20011 ) considered a simple system in which H and O atoms land on classical 



grains to produce H2, OH, and O2. Biham and co-worke rs focused on the formation of H2 



With a slow diffusion rate for H based on t he analvsis bv 



Katz et al 



work by Pirronello, Vidali, & co-workers (jPirronello et al 



( 19991 ) of experimental 



1997al )br). they found the mas- 



ter equation approach to be needed mainly for small grains. For the formation of H2, the 
dynamical variables are the probabilities Ph(A^h) that there are Ah atoms on the grain at 
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time t. The time derivatives Ph(^h)) = 0, 1,2, ... , are expressed in terms of the ad- 
sorption, reaction and desorption terms. The master equation provides the time evolution of 
Ph(^h), = 0, 1, 2, . . ., from which the recombination rate can be calculated. The master 



equation has also been used for the study o: 



fStantcheva et al. 



2002; 



• react ion networks involving multiple species 



Stantcheva fc Herbst 



20031) . 



In addition to th e master equatio n method, a Monte Carlo method (|Tielens fc Hagen 



19821 1 was propo sed (ICharnleyl I2OOII), and applied to surface reaction networks involving 



multiple species fjCaselli et al 



3 



200^ Stantcheva fc HerbstI 



20031 ). The weakness of the Monte 



Carlo approach to surface chemistry is that it cannot be used easily for models in which gas- 
phase reactions take place simultaneou sly. Finally, a semi- empirical approach, known as the 



modified rate method, was proposed by 



Caselli et al 



( 19981 ) . Althoug h artificial, the modifiec 
rate method is easy to employ and has been used with mixed success f Stantcheva et ahlboOl 



Caselli et al. 



20021) . 



In addition to H, which is the most abundant atom in the universe, there is a density of 
D atoms, which is smaller by a factor of about 10"^ compared with the density of H atoms 
in regions where the gas is primarily atomic. If deuterated hydrogen molecules, namely HD 
and D2, are formed on grain surfaces in the same way as H2, these molecules are expected 
to be rare due to the low abundance of D atoms. Of course, some differences between the 
formation rates of HD and D2 and that of H2 are expected because of the possibly different 
diffusion and desorption rates of H and D atoms on the surface. It was realized early in 
the study of interstellar chemistry, however, that the formation of HD in diffuse interstellar 
clouds occurs more efficiently in the gas phase via the sequence of reactions 



H" 
D- 



+ D- 
+ H2 



D+ + H 
^ HD + H+, 



(1) 



1973! ) and explain the observed HD/H2 abundance ratio in diffuse clouds. By the stage in 



which can lead to enhancements of ~ 10^ over surface formation ()Watsonll973 ; 



Dalgarno et al 



which dense interstellar clouds have been formed, HD has become the major repository 
of deuterium, and is the starting point for a rich che mistry of deuterium fractionation 



that occurs in the g a s phase via ion-rnolecule r eactions (IWatson 



1997 



Roberts fc Millar 



2000aOb 



I 



Rodgers &: Charnlev 



Millar et a1 



2001 



IM 



Roberts et al 



Charnlev et al 

2002! ). Since gas-phase fractionation produces a large D/H atomic abundance ratio in dense 
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clouds, it then becomes possible that surface processes can produce multiply deuterated 



species in high abundance following the adsorption of D atoms ( Charnlev et alll997l : modgers fc Charnlev 



Caselli et al 



2mA 



Sta,ntcheva Herbs 



3|2q3). 



Although it has been thought for some time that the formation of HD in diffuse clouds 
occurs primarily in the gas, the surface processes for the formation of HD and D2 have 
not been looked at carefully. In this paper we propose a mechanism taking place on the 
surfaces of dust grains that, under certain conditions, can provide a significantly enhanced 
production of deuterated molecules. We demonstrate the mechanism for the simple case 
that involves only H and D atoms in diffuse clouds, and examine the formation rates of HD 
and D2 vs. that of H2. The mechanism is based on the assumption that D atoms stick to 
the surface more strongly than H atoms so that their desorption rate is lower. An isotope 



effect of this type ha s been observed in various experimental situations (jKoehler et al, 



Hoogers et al. 



mm- 



Due to the very low cosmic abundance ratio of D vs. H (10~^) most of the atoms that 
hit the grain surfaces in diffuse clouds are H atoms. However, the residence time of D atoms 
on grains is expected to be longer than that of H atoms, due to their lower desorption rate. 
Thus, the D/H atomic abundance ratio on the grains is enhanced compared with the gas- 
phase ratio. As a result, newly adsorbed H (or D) atoms are more likely to find D atoms 
already residing on the grains, giving rise to an enhanced production of HD molecules. This 
enhancement turns out to be even larger on very small grains on which only a handful of H 
and D atoms typically reside. This effect is expected to be significant in diffuse clouds, where 

the grain size distribution is broad and dominated 



3 phase is lov 


/ anc 


Li & Draine 


2001) 



The paper is organized as follows. The interaction between hydrogen atoms and dust 
grains in interstellar clouds in described in Sec. 2. The formation of molecular hydrogen on 
large grains is described in Sec. 3 using the rate equation model. The formation of molecules 
on small grains is presented in Sec. 4 using the master equation approach, followed by an 
asymptotic analysis for very small grains in Sec. 5. The simulations and results are presented 
in Sec. 6. The results are discussed in Sec. 7 and summarised in Sec. 8. In Appendix A we 
consider the relevance of the isotope effect studied here to related experiments. 
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2 THE INTERACTION BETWEEN HYDROGEN ATOMS AND DUST 
GRAINS 

Consider a diffuse interstellar cloud dominated by a density (cm~^) of H atoms and also 
including a density ud = fuu of D atoms, where / is the D/H abundance ratio (the cosmic 
ratio is approximately / = 10^^). The typical velocity of the hydrogen isotope I (H or D) in 
the gas phase is denoted by vi (cm s~^). It is given by 



(2) 



/ 8 fc^Tgas 

Vl = \ —, 

V vr mi 

where the mass mi is mu = 1.67 x 10~^^ (grams) for H atoms and mo = 3.34 x 10"^"^ (grams) 
for D atoms and the temperature of the gas is denoted by Tgas. The cloud includes a density 
ngr (cm^^) of dust grains. To evaluate the flux of atoms onto grain surfaces we will assume, 
for simplicity, that the grains are spherical with a radius r (cm). The cross section of a grain 
is 0" = vrr^ and its surface area is Aa. The number of adsorption sites on a grain is denoted 
by S, and their density s (sites cm~^) on the surface is given by s = S/{Aa). The fluxes 
Fu (atoms s~^) and Fd of H and D atoms onto the surface of a single grain are given by 
Fi = uiVia. The atoms stick to the surface and hop as random walkers between adjacent 
sites until they either desorb as atoms or recombine into molecules. The desorption rates of 
H and D atoms on the surface are given by 



Wi = u-exp[-E,{l)/kBT], 



(3) 



where z/ is the attempt rate (typically assumed to be 10^^ s~^), £"1(1) is the energy needed 
for desorption of an atom of isotope I and T is the surface temperature. The hopping rates 
of the atoms from site to site are 



ai = u- exp[-Eo{T)/kBT], 



(4) 



where Eq{T) is the activation energy barrier for diffusion of the isotope I. Here we assume 
that diffusion occ urs only by thermal hopping, in agreement with recent experimental results 



( Katz et al. 



1999) 



allowed to occur (jHasegawa et al 



. In the faster diffusio n rates often adopted by astrochemists, tunneling is 



19921 ) and results in very different rates for H and D. The 



rate Ai = ai/S is approximately the inverse of the time ts{l) required for an atom of isotope 
I to visit nearly all the adsorption sites on the grain surface, because in two dimensions the 
number of distinct sites visited by a ra ndom walker is linearly p roportional to the number 
of steps, up to a logarithmic correction (|Montroll fc Weisa 119651 ). 
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In this p aper, we take the density of adsorption sites on the surface as s = 5 x 10^^ 
(sites cm"^) (jBiham et aLll2nnih . The density of H atoms is riH = 50 (cm~^) and the density 
of D atoms is ud = frin, where / = 10~^ is the cosmic abundance ratio, which according 
to Eq. (0) gives rise to to the ratio Fd/Fh = 10^^ /V2. The gas temperature is taken as 
Tgas = 90 K, thus vh = 1-37 x 10^ and vd = 9.7 x 10^ (cm s~^). The surface temperature is 
T = 18 K, although variations from this value are considered. These values are within the 
range of typical parameters for diffuse interstellar clouds. 

The parameters used for the interaction between the adsorbed H and D atoms and the 
grain are crucial to the determination of the isotope effect. Since D is heavier than H, it sits 
lower than H in any potential surface when zero-point energy is included. The smaller zero- 
point energy of D leads to larger values for the desorption energy and the barrier against 
diffusion, although the differences between H and D are uncertain. In the analysis of their 



experimental results, for example, iKatz et al 



D atoms at all. In a previous paper. 



(Il99fl) assu med no difference between H and 



Caselli et al 



( 20021 ) argued for a difference of 2 meV 



between both the desorption energies and barriers against diffusion. Here, the adopted energy 
barriers for diffusion are Eq{H) = 35 meV for H atoms and Eq{D) = 37 meV for D atoms, 
while the desorption energies are Ei{H) = 50 meV for H atoms and Ei{D) = 60 meV for 
D atoms. The values for H atoms are intern iediate between t hose obtained experimentally 
on olivine and amorphous carbon surfaces ( Katz et al. 19991 ). The difference between the 



desorption energies is assumed to be 10 meV; an argument in favour of this larger value is 
that the barrier against diffusion balances the zero-point energy of a potential well with that 
of a saddle point while desorption does not possess a saddle point. In any case, the difference 
in desorption energies between H and D is crucial to the isotope effect discussed here. 



3 MOLECULAR HYDROGEN FORMATION ON LARGE GRAINS 

After being adsorbed onto the surface, the H and D atoms hop between adsorption sites 
until they either desorb or meet one another and form molecules. The number of atoms of 
isotope I on the surface of a grain is denoted by A^i. The rate equation analysis deals with the 
expectation values (A^i), I = H, D obtained by averaging over a large ensemble of identical 
grains under given physical conditions. The rate equations for the expectation values of the 
H and D populations take the form 
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d{Nu) 



Fu - W^{Nn) - 2An{Nuf - (Ah + Aj,) (Nu) (Nr,) 



dt 



d{Nj,) 



Fd - W^D(iVD) - 2A^{Nj,f - {An + A^){Nn){N^). 



(5) 



dt 



These equations provide the time derivatives of the sizes of the populations of reactive 
species. Each equation includes a flux term, a desorption term and two reaction terms. 
The first reaction term accounts for reactions between atoms of the same specie, while the 
second reaction term accounts for the reaction that produces HD molecules. The rate of the 
latter reaction is proportional to the rate at which H and D atoms encounter each other, 
which is proportional to the sum of the hopping rates of the two species. Therefore, the 
pre-factor of this term is the sum of An and Ad. By integration of the rate equations for 
a given initial condition one obtains the time dependence of the population sizes. The rate 
equations provide a good description of the surface reactions with our chosen parameters as 
long as the grain is not too small and the flux is not too low. The rate equations are strictly 
valid when both the D and H populations per grain are greater than unity. However, we 
have found that, despite a very small D population in even our largest grains, as long as 
the H population per grain is not too much smaller than unity, the rate equations provide 
a reasonable approximation to the master equation results. 

We will now solve the rate equations for the case of steady state conditions, in which 
d{Ni)/dt = for I = H, D. Due to the small ratio between the fluxes of D and H atoms, 
the effect of the D atoms on the density (Nn) can be neglected. Thus, neglecting the HD 
formation term in the first equation in Q we obtain 




(6) 



Inserting this solution into the second equation in Q we obtain 




(7) 



where 



W,s = WD + {Au + AD){Nn) (8) 

is the effective desorption coefficient for the D atoms. The ratio (A^d)/(^h) between the 
population sizes of D and H on the grain surface is obtained from Eqs. (jH)) and ((Tj). This 
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ratio turns out to be higher than the gas phase ratio D/H due to the lower desorption rate 
of D atoms. The production rate of HD molecules per grain is 

i?HD = (^H + ^D)(A^H)(A^D) (9) 

and of H2 molecules is 

i?H, =^H(iVH)^ (10) 

The ratio between these production rates is thus 

Rhd _ ^1 + (M. (11) 



Due to the assumption that Eq{H) < Eq{D), namely that the diffusion of H atoms is faster 
than of D atoms, we obtain that the additional enhancement factor multiplying (A^d) / (-^h) 
is in the range 1 < 1 + An/A-a < 2. The production rate of D2 molecules is 

Rb, - Aj,{Nj,)^, (12) 



(13) 



and the ratio between the D2 and H2 production rates is 
Rd,^A^{N^ 
Rn, ^H(iVH)2- 

4 MOLECULAR HYDROGEN FORMATION ON SMALL GRAINS 

The master equation that describes the hydrogen-deuterium system consists of a two- 
dimensional matrix of equations. These equations describe the time derivatives of the prob- 
abilities -P(A^H, -^d) of having A^h hydrogen atoms and A^d deuterium atoms simultaneously 
on the surface. The master equation takes the form 

P(Nu,Nj,) = Fn[P{Nji-l,NT,)-P(Nji,Nj,)] 

+ F^[P{Nn,Nj,-l)-P{Nn,Nu)] 

+ W^h[(A^h + 1)P(A^h + 1,A^d)-A^hP(A^h,A^d)] 

+ W^d[(A^d + 1)P(A^h,A^d + 1)-A^dP(A^h,A^d)] (14) 

+ An[{Nn + 2){Nn + l)P{Nn + 2, Nn) - Nu{Nn - l)P{Nn, Nq)] 

+ Ag[{N-d + 2){Nj, + 1)P(A^H, A^D + 2) - Nj,{N-D - l)P{Nti, Nj,)] 

+ (An + Ad)[(A^h + 1)(A^D + 1)P(A^H + 1, A^D + 1) - NjiNr^PiNn, A^d)], 

where A^h, A^d — 0, 1, 2, Each of these equations includes three sets of terms. The first set 



Enhanced formation of HD and D2 on dust 9 
describes the effect of the incoming ffuxes Fh and on the probabihties. The probabihty 
P(A^H, ^d) increases when an H atom sticks to the surface of a grain that already has A^h — 1 
adsorbed H atoms and A^d adsorbed D atoms. There is a similar term for the adsorption of a 
D atom. The probability P{Nii, Nj)) decreases when a new atom is adsorbed on a grain with 
A^H H atoms and A^d D atoms that are already adsorbed. The second set of terms describes 
the effect of desorption. An H [D] atom that is desorbed from a grain on which there are 
A^H adsorbed H atoms and A'^d adsorbed D atoms decreases the probability P(A^H; ^d) and 
increases the probability P(A^h — 1, A^d) [P(A^h, A^d — !)]• The third set includes three terms 
describing the reactions that lead to the formation of H2, D2 and HD, respectively. The 
rate of each of these reactions is proportional to the sum of the sweeping rates of the atoms 
involved (Ah or Aq) times the number of pairs of atoms that may participate in the reaction. 
In the equations in which A^h = 0, 1 or Nj) = 0, 1, some of the terms vanish. These terms 
can be easily identified since their initial or final states include a negative number of H or 
D atoms on the grain. 

By suitable summation over the probabilities, the master equation provides the moments 

of the population size distribution, given by 

00 00 

(A^hO= E E «P(iVH,ArD), (15) 

where k,i = 0,1,2,... [note that in case that k = {i = 0) the corresponding power, 
(^d); is removed from the equation]. The order of each moment is specified by the 
sum of these powers, namely hj m = k + i. Clearly, there are two first-order moments, 
(Ah) and (A^d), three second-order moments, (A^h); (^d) (^h^d), as well as m -|- 1 
moments of any higher order m > 2. The average rates of production of the H2, HD and 
D2 molecules on the surface of a single grain can be obtained by a proper summation over 
the third set of terms in the master equation (fT^ . These production rates can be expressed 



in terms of the first and second m oments of the population size distribution (jBiham et al 



2001 



Stantcheva et al. 



2001 



20021 ): 



Rn, = An{{Nl)-{Nu)) 
Rkd = {An + Aj,){NnNu) 

Rj,, = Ad ((Ar2) - (ATd)) . (16) 
The fact that the production rates depend on the second moments explains the failure of the 
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rate equations in the limit in which the average surface populations of reactive species become 
^ 1. In the rate equations, which account only for the first moments, second moments are 
approximated by products of the corresponding first moments. The equations fail to account 
for the production rates when this approximation is no longer valid. 

The production rates given by Eqs. ()16|) can be obtained by direct integration of the 
master equation (fT^ as a set of time dependent coupled differential equations, using a 
standard Runge-Kutta stepper. In this paper we focus on the production rates under steady 
state conditions. Therefore, the data collection is done only after the populations of the 
reactive species have reached their asymptotic sizes under the simulated conditions. In order 
to limit the number of equations, suitable cutoffs are imposed, either on the population 
size of e ach reactive specie or on the entire population of adsorbed atoms on the grain 



surface (Stantcheva et al. 



20011 ). The production rates can also be obtained from Monte 



Carlo simulations, based on the master equation (jCharnlev 



mm- 



5 THE LIMIT OF VERY SMALL GRAINS 

In the limit of very small grains, one can reach the situation in which both A^h ^ 1 and 
A^D -C 1. In this limit the formation rate of molecules is reduced and the dominant mechanism 
for removing atoms from the surface is thermal desorption. The formation of molecules can 
then be considered as a perturbation on top of the adsorption and desorption processes. 
This limit can be described by the master equation using cutoffs that do not allow more 
than a single atom on the grain. The master equation then consists of only three coupled 
equations, for P(0,0), -P(1,0) and P{0, 1). With these cutoffs, the average population sizes 
are related to the probabilities by 



(iVn) = P(1,0) 

(iVc) = P(0,1). (17) 
Solving the three coupled equations for steady state conditions we obtain that 



(Nn) = Fu/Wu 

Consider the situation in which there is one atom on the surface of a grain. 



(18) 

Since the grain 
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is very small an atom can easily scan the whole surface before it desorbs, namely Ai^Wi. 
Thus, when a second atom is adsorbed, it is most likely to form a molecule with the atom that 
is already on the surfac e. Therefore, the produc tion rate does not depend on the diffusion 
rates, a point noted by iTielens fc Haeenl (|l982^ and subsequently used to define the term 
"accretion" limit rigorously, as distinguished from the "reaction" limit. In the former limit, 
the production rates are given by the flux terms of H and D atoms in the master equation, 
multiplied by the probabilities of having a single atom (either H or D) on the surface, namely, 
FiP(0, 1) and -FiP(l, 0), where I = H, D. Using Eqs. (fTTj) . the production rates of molecules 
take the form 



Rub = {Nn)F^ + {Nj,)Fn 

Rb, = {No)Fo. (19) 

The HD molecule can form in two channels: an H atom on the surface reacts with a newly 
adsorbed D atom, or a D atom on the surface reacts with a newly adsorbed H atom. Due 
to the longer residence time of a D atom on the surface, the second channel is dominant. 
Therefore, the production ratio Rud/Rr2 is approximately equal to the ratio between the 
population sizes (A^d)/(^h) on the surface. This amounts to an enhancement of the HD 
production by a factor of Wu/Wd compared to the flux ratio, namely 

-^HD ^ WnFh /2Q\ 
Rh, WnF^- ^ ^ 

The production rate of HD is thus enhanced with respect to that of H2 by a factor given 
by the ratio of the desorption rates of H and D atoms. Since the desorption rates depend 
exponentially on the activation energy barriers for desorption Ei{H) and Ei{D), a small 
difference between these barriers can produce a large enhancement in HD production. The 
enhancement in D2 production is given by 



Rb, Wn (Fj, 



Rr, Wj, VFh 



(21) 



giving rise to a similar enhancement compared to the expected production rate based on the 
ratio between the fluxes of D and H. 
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6 SIMULATIONS AND RESULTS 

To evaluate the production rates of H2, HD and D2 molecules we have used the master 
equation approach given by Eq. ()14j) . The results presented in this paper were obtained by 
direct numerical integration of the set of differential equations until they reached steady 
state conditions. For a range of very small grains, namely 5* < 5000, we also performed 
kinetic Monte Carlo simulations based on the master equation, and found perfect agreement 
between the two methods. 

In Fig. 1 we show the average population sizes of the atomic species on the grain surface, 
(A^h) (solid line) and (A^d) (dashed line), as a function of the number of adsorption sites, 
S, and grain radius, r, under steady state conditions. The results are presented for a range 
of small grains on which the average number of H atoms is smaller than one while the 
average number of D atoms is of the order of 10~^ and less. Under these conditions the 
master equation is required in order to evaluate the production rates correctly, although 
the rate equations (dotted lines) are progressively better as the size of the grain increases. 
Specifically, one can see that the rate equation approach appears to handle the average 
hydrogen abundance well for all grain sizes, while becoming much worse for the average 
deuterium abundance at small grain sizes. One should note that the actual production rate 
of molecules does not depend on (A^h) and (A^d) only [see Eq. (fT6|l ]. so one cannot argue 
about the suitability of the rate equation method for the determination of production rates 
based on the results shown in Fig. 1. 

The production rates of H2, HD and D2, molecules per grain vs. 5* and r are shown in 
Fig. 2. Clearly, the production rates decrease as S is reduced. For large grains that support 
sizable populations of H and D, and can be described by rate equations, the production rates 
per grain are proportional to the surface area, namely linear in S. In the limit of diminishing 
grain sizes, the production rates become quadratic in S, as can be seen by substitution of 
(A^h) and (A^d) into Eq. (fT9|) . which leads to a quadratic dependence on flux. The range of 
grain sizes presented in the simulation results reaches both the limit of large grains and the 
limit of small grains, covering the intermediate range between these two limits. 

In order to quantify the relative enhancement in the formation rate of deuterated molecules 
vs. H2 we present in Fig. 3 the ratios -Rhd/-Rh2 and Rd^/ R^^ vs. S and r. It is observed that 
as S decreases these production ratios increase by over a factor of 20. In the limit of large 
grains, the production ratio converges towards the results of the rate equations (dotted hne 
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on right-hand side). In the hmit of very small grains it converges towards the results of the 
asymptotic analysis presented in Sec. 5 (dotted line on left-hand side). 

To examine the dependence of the efficiency of the molecular formation on the grain size 
we show in Fig. 4 the production rates per site of H2, HD and D2 on the grain as a function 
of 5* and r. In the limit of large grains this quantity saturates and does not depend on the 
grain size. However, as the grain size decreases, the production rate per site goes down. The 
production rate of H2 is reduced more quickly than that of the deuterated species. This gives 
rise to an enhancement in the relative production rate of HD and D2 vs. H2. 

To account for the production rate of molecular hydrogen in i nterstellar clouds one should 



integ r ate oyer the entire distribution of grain sizes in the c loud (jMathis et al 



1990 



1996 



Draine fc Lee 



1984; 



Draine &: Anderson 



1977 



Mathis 



19851). Consider a cloud in which the 
grain size distribution is given by ^gr(^) (cm^"^) in the range r^nin < r < rmax, where rmin and 



max! 

~3^ 



'"max are the lower and upper cutoffs of the distribution [here ng^{r)dr (cm ) is the number 
density of grains of sizes in the range (r,r+dr)]. The formation rates TZr^, 7?.hd and 7^d2) 



(cm s ^ ) in the cloud are given by 



and 



ngr{r)Rji^{r)dr, 
ngr.{r)R^Y){r)dr, 

ngr{r)RD2ir)dr, 



(22) 
(23) 

(24) 



respectively, where Ru_2{r), RnDi^) and Rd2{^) ^-re the production rates per grain for grains 
of radius r. 

O bservations support the assumption tha t the grain size distribution takes the power-law 



form ( Mathis et al 



ma 



Weingartnei 



2(M) 



in a range of sizes r^^ 

(a - l)ngr(total) 



where 



K 



min max 



(25) 



(26) 



and ngr(totai) (cm~^) is the total density of grains, integrated over the whole range of grain 
sizes. To exemplify the role of the grain size distribution we analyze the case in which 
a = 3. This distribution has the special property that the mass of grains in the cloud is 
equally distributed within the range of grain sizes. The total surface area of all grains is 
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distributed such that grains in the size range between (r, r + dr) have a total surface area 
proportional to 1/r. As a result, most of the surface area belongs to grains in the lower end 
of the size distribution. In the analysis we chose rmin = 4 x 10"'' cm and rmax = 2 x 10~^ 
cm. The parameter K is chosen such that the density of large grains, of radii in the range 
5 X 10^^ < r < 2 X 10~^ cm, is (classical) = lO^^^nn. Integrating the grain size distribution 
over the entire range we obtain ngi.(totai) = 1-66 x lO^^'^nn- The minimum grain size is 
chosen to be larger than PAH molecules so that the problem of stochastic heating and 
desorption bv photons doe s not domi nate normal therra al evaporation. According to the 
calculations of lDraine fc Lil ()2nnih and lLi fc Draind 1)20011 1. however, grains of our minimum 
size may experience stochastic heating of temperatures up to 40 K at time intervals (~ 10^ 
s~^) somewhat shorter than the D evaporation time at the grain temperature chosen. Such 
heating would effectively shorten the grain lifetime of D atoms and reduce the isotope effect. 
The time scale between heating events could be lengthened by extinction in the cloud. 

To examine the contribution of grains of different sizes to the total production rates of 
H2, HD and D2 we show in Fig. 5 the differential production rates given by the integrands 
of Eqs. and as a function of the grain radius r. Two competing effects play a 
role in this process. The production efficiency of molecules on the grain decreases as the 
grain size is decreased. However, small grains account for most of the surface area and thus 
provide a large fraction of the molecules that are produced. The relative enhancement in the 
production of deuterated species arises because their production efficiency decreases more 
slowly than that of H2 as the grain size is reduced. 

Integrating over the range of grain sizes we obtain that the production rates are 7^h2 = 
3.9 X 10^^^ cm-3 s-\ 7^HD = 1-2 x 10"^^ cm^^ and 7^D2 = 1.5 x lO^^s ^^^-3 g-i_ 
Therefore, the ratio 7?.hd/'^H2 = 3.1 x 10^"^ (71d^/71ii2 = 3.9 x 10^^) is larger by about 
a factor of 44.1 (78.3) than the expected value based on the ratio of the incoming fluxes, 
namelv Fn/F^ = 1 0~^/a/2. Observations indicate that the exponent in Eq. ()25|) is a = 3.5 



( Weingartner 



20011 ). Therefore, the contribution of very small grains may be even more 



dominant, leading to a larger enhancement in the formation of HD and D2 vs. H2 than in 
the analysis shown here. 

We will now compare our results with the standard rate law used for the production 
of molecular hydrogen in diffuse regions. This rate law assumes that all the grains are of 
radius r = 10~^ cm and that their density is equal to ngi.(ciassicai)- It is also assumed that the 
efficiency is rj = 1, namely, all the atoms that collide with grains end up in molecules. Under 
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these assu mptions the production rate of H2 molecules is given by T^-Ha = ''^H'yH^gr(ciassicai)C"/2 



cm ^ s ^ ( Hollenbach et al 



19711 ). or more conveniently, by 



7^H2 = R{^2) nn n, (27) 

where rin is the density of H atoms and n is the total density of H nuclei in all forms 
(although we assume here that all the hydrogen in the gas is in the atomic form, namely 
n = nn). The parameter -R(H2) (cm^ s~^) is the effective rate coefficient given by 

i?(H2) = 5 X lO-^yTgas/300. (28) 

For the conditions considered here, -R(H2) = 2.74 x 10^^^ cm^ s~^ and TZr^ — 6-85 x 10~^^ 
cm~^ s^^. This result is about an order of magnitude higher than our result, which is 
consistent with R{B.2) = 1.56 x 10^^^ cm^ s~^. The reason for this apparent discrepancy 
is the low production efficiency of H2 molecules under these conditions, which on large 
grains is rj = 0.09 (given by the rate equations) and is even lower for small grains. The 
production rate of HD molecules (cm~^ s~^) can be expressed in a similar way, by 

7^HD = i?(HD) no n, (29) 

where i?(HD) (cm^ s~^) is the effective rate coefficient for this reaction. The values of the 
rate coefficients -R(H2) and i?(HD) under the conditions studied here, are shown in Table 1, 
for a = 3 and 3.5 and for several grain temperatures between 12 and 20 K. 

In Fig. 6 we present the temperature dependence of the production rates of H2 and 
HD molecules. The production rates TZr^: obtained from the rate equations, are shown 
by the solid lines for the case that the exponent of the grain size distribution is a = 3.5 
(upper solid line) and for a = 3 (lower solid line). The master equation results for these 
two cases are also shown (by o and x, respectively). The production rates T^hd, obtained 
from the rate equations, are shown by the dashed lines, for a = 3.5 (upper dashed line) 
and for a = 3 (lower dashed line). The master equation results are also shown (by o and 
X, respectively). These graphs provide much insight on the enhancement of HD production. 
They show that for grain temperatures lower than about 16 K, the production efficiency is 
7] ~ 1 and the rate coefficient -R(H2) is consistent with Eq. ()28|) . At this temperature range 
there is no enhancement in the production of HD molecules. The enhancement appears in 
a range of temperatures, roughly between 16 K and 20 K, where the H2 production rate 
quickly decreases as T increases while the HD production rate remains steady. Above 20 K 
the production rates of both H2 and HD are diminished. 
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In Fig. 7 we present the ratio, 7^hd/'^H2) of the production rates of HD and H2 molecules 
as a function of the grain temperature, as obtained from the rate equations (solid line), and 
the master equation for a = 3 (x) and a = 3.5 (o). The enhancement ratio increases sharply 
as the temperature is raised above 16 K and reaches a peak around 20 K. A significant 
enhancement is observed on large grains, for which the rate equation results apply. When 
the grain size distribution is taken into account, further enhancement is obtained due to the 
contribution of the small grains, on which the enhancement ratio is even larger. 

The dependence of the ratio of the production rates -Rhd/-Rh2 on the grain size 5* and 
radius r is shown in Fig. 8 for several grain temperatures. At 14 K there is no significant 
enhancement for any grain size. At 16 and 18 K the enhancement strongly depends on 
the grain size and increases sharply as the grain size decreases. At 20 K there is a large 
enhancement for grains of all sizes. However, at this temperature the total production of 
both H2 and HD is very small, thus making little contribution to the gas phase abundances. 



7 DISCUSSION 



In this paper we have used a simple description of the grains as spherical objects. In practice 
the grains are expected to be amorphous with rough surfaces and may include pores. These 
features may increase the surface area of the grain and introduce a distribution of diffusion 
and desorption barriers. However, no solid data are available to quantify such effects. More- 
over, the analysis of experi mental results on molecu la r hydr ogen production on olivine and 



amorphous carbon surfaces (IPirronello et c 



energy barrier for each process (jKatz et al 



1997ai bl. Il999[ ) did not require more than one 



The analysis presented in this paper is restricted to the case in which the mobility of 
adsorbed atoms on the surface is due to thermal hopping. Tunneling of H and D atoms 
on the surface is not included since experimental results indicate that tunneling does not 
play a signi ficant role in the format io n of m olecular hydrogen under astrophysically relevant 



conditions ( Pirronello et al 



1997al jbl. I1999I ) . In a recent anal ysis of the results of the se ex- 



periments, chemisorption processes were taken into account (ICazaux &: TielensI 120021 ). The 
study of chemisorption processes is motivated by observations of molecular hydrogen in 
photon-dominated regions in which the typical grain temperature is around 40 K, making 
the production of H2 molecules by physisorbed H atoms hopelessly inefficient. 

The mechanism presented here for an enhanced production of HD and D2 molecules, vs. 
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H2, takes place in a window of grain temperatures. For the surface parameters used here 
this temperature window is between 16 - 20 K. In this temperature range the production 
efficiency of H2 is reduced, while the efficiency of HD and D2 production remains steady. 
The effect is observed for grains of all sizes, but the relative enhancement is larger on small 
grains. It will be interesting to see if the enhancement plays a significant role in chemical 
models of diffuse and translucent clouds; relevant model calculations are planned by Le Petit 
and Roueff (private communication). 



8 SUMMARY 



The formation of H2, HD and D2 molecules on dust grain surfaces was studied using the 
master equation approach. The production of these molecules was assumed to occur via the 
Langmuir-Hinshelwood (diffusive) mechanism. It was shown that for surface temperatures 
in excess of 15 K, the fractions of HD and D2 molecules in the population of hydrogen 
molecules that are produced are significantly enhanced compared with their expected ratio 
based on the gas phase abundances. The enhancement takes place on grains of all sizes, but 
is largest on the smallest grains considered until the temperature reaches 20 K. It arises 
because of the assumption that D atoms are bound slightly more strongly to grains than 
are H atoms. At a surface temperature of 18 K, with binding energies and rates of diffusion 
based on laboratory experiments, the enhancement is found to be a factor of 44.1 (78.3) 
for the production of HD (D2)with the assumption that the grain size distribution goes as 
r~^. Although the enhancement is even larger at higher temperatures, the absolute rates of 
production for H2 and its isotopomers become too small to be of importance. Use of the 
master equation method can lead to somewhat greater enhancements (factor of a few) than 
the simple rate equation method.^ 



While gas-phase mechanisms (jWatson 



19731 ) may account for the observed abundance of 



HD, it is possible that surface mechanisms such as the one presented here also play a role 
in the production of deuterated molecules. Recent observations have uncovered a variety of 
multiply deuterated inolecular species in the in terstellar medium such as D2CO and ND3 



( Ceccarelh et al, 



mm 



van der Tak et al, 



20021 ). Gas phase reaction networks that lead to 



the formation of such mu 



examined ( Charnlev et al 



tiply deuterated species have been studied and the i r effic iency 



1997 



Rodgers fc Charnlev 



2001 



Roberts et al. 



2002 



20031) . The 



picture in pre-stellar cores is that gas-phase fractionation can explain much of what is 
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observed, but the picture in protostellar cores, where grain mantles tend to be evaporated 
back into the gas, is that surface r eactions plav an impor tant role in the fractionation 



leading to multi-deuterated species (jStantcheva et al 



20021 ) ■ It will be interesting to see 



if the surface mechanism presented here may be extended to pro duce multiply deuterated 



species. Preliminary work o n this matter has been undertaken by 



Caselli et al 



(2002) and 



Stantcheva fc HerbstI (|2003h . 
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Table 1. The rate coefficients R(H2) and i?(HD) (cm^ s ^) for several grain temperatures between 12 and 20 K and for grain 
size distributions characterized by a = 3 and 3.5. 

T R{B.2) R(HD) 



a = 3 a = 3.5 — 3 cx = 3.5 

12 4.31(-17) 8.35(-17) 6.1(-17) 1.18(-16) 

14 4.40(-17) 8.18(-17) 6.36(-17) 1.2(-16) 

16 2.29(-17) 3.03(-17) 6.33(-17) 1.19{-16) 

18 1.56{-18) 1.67(-18) 4.85(-17) 7.53(-17) 

19 2.22(-19) 2.46(-19) 2.82(-17) 3.59(-17) 

20 3.11(-20) 3.92(-20) 7.19(-18) 9.19(-18) 
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Figure 1. The average number of H atoms, (A^h) (solid line), and the average number of D atoms, (A^d) (dashed line) on the 
surface of a grain that has S adsorption sites, as a function of S and grain radius r, obtained from direct integration of the 
master equation. The corresponding rate equation results are shown by dotted lines. 




Figure 2. The production rates per grain of H2 (solid line), HD (dashed line) and D2 (dashed-dotted line) molecules as a 
function of the number of adsorption sites, S, on the grain surface and grain radius r. In the limit of large grains the production 
rates approach a linear dependence on 5, while in the limit of very small grains they approach a quadratic dependence on S. 
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Figure 3. The ratio between the production rates of HD [D2] and H2 molecules, Rhd/ Rh-^ (dashed line) [Rd^/ Rh^ (dashed- 
dotted line)] vs. the number of adsorption sites S on the grain surface and grain radius r. In the limit of large grains both ratios 
converge to the rate equation results (dotted lines on the right-hand side). As the grain size goes down, these ratios increase 
by over a factor of 20 until they converge to the results of the asymptotic analysis for the limit of very small grains (dotted 
lines on the left-hand side). 
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Figure 4. The production rates of H2, HD and D2, molecules per site on the surface of a grain that has 5 adsorption sites, as 
a function of S and grain radius r. These curves quantify the dependence of the efficiency of molecular production per surface 
area on the size of the grain. While the production rate of H2 declines sharply as the grain size decreases, there are more 
moderate reductions in the production rates of HD and D2 . 
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Figure 5. The differential production rates of H2, HD and D2 molecules [given by the integrands in Eqs. I22i . I23i and 1241 . 

respectively] on the surfaces of grains of radii in the range r, r + dr. 
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Figure 6. The production rates, obtained from the rate equations, of H2 (sohd Unes) and HD (dashed lines) molecules as a 
function of the grain temperature, for a = 3.5 (upper line in each pair) and for a = 3 (lower line in each pair). The master 
equation results for these two values of a are also shown (by o and X , respectively) . 
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Figure 7. The ratio between the production rates of HD and H2 molecules as a function of the grain temperature, as obtained 
from the rate equations (solid line), and the master equation for a = 3 (x). and a = 3.5 (o). The enhancement ratio increases 
sharply as the temperature is raised above 16 K and reaches a peak around 20 K. The enhancement exists for all grain sizes 
with further enhancement in the case of small grains. 
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Figure 8. The ratio of the production rates, i?HD/KH2i grain size S and radius r, for several grain temperatures. At 
T = 14 K there is no significant enhancement. At T = 16 and 18 K the enhancement strongly depends on the grain size and 
increases sharply as the grain size is reduced. At T = 20 there is a large enhancement for grains of all sizes. However, at this 
temperature the total production of both H2 and HD is very small. 



Enhanced formation of HD and D2 on dust 31 
APPENDIX A: ISOTOPE EFFECTS IN RELATED EXPERIMENTS 



Isotope effects resembling the type studied he re were obse rved experimentally on various 
surfaces including metals (iHoogers et al.iiil995i) and silicon (|Koehler et al.l 119881 ). However, 
these experiments were done at much higher temperatures than those relevant for molecular 
hydrogen formation. Also, the surfaces on which isotope effects were studied are more or- 
dered than the astrophysically relevant surfaces. Therefore, it may not be possible to draw 
conclusions from these experiments for astrophysically relevant surface processes. 

A crucial feature in the recent experiments on molecular hydrogen formation is the 
use of two beams, of H and D atoms, and the measurement of the HD production rate 



( Pirronello et al 



1997a b 



1999! ). The use of HD molecules improves the signal-to- noise ratio 



in the experiment and allows reliable measurement of the production rate. In the analysis of 
the experimental results it was assumed, for simplicity, that there is no isot ope effect, iiamely , 
that the diffusion and desorption barriers for H and D atoms are identical ( Katz et al.lll999 ). 
This assumption was made in order to limit the number of free parameters used in the 
fitting. Under the experimental conditions this assumption had little effect on the results of 
the analysis. To explain this issue, let us recall the experimental procedure. The experiment 
followed the standard temperature programmed desorption (TPD) procedure, which consists 
of two stages. In the first stage, the surface was irradiated by two beams: one of H atoms 
and the other of D atoms. The beams were of (approximately) equal fluxes and lasted for 
an equal irradiation time, giving rise to equal coverages of H and D atoms. In the second 
stage the surface temperature was raised at a constant rate and the production rate of HD 
molecules was recorded using a mass spectrometer. 

The analysis of the experimental results was based on the fact that the surface was of 
macroscopic size, so that rate equations of the form (jS)) can account correctly for the atomic 
populations on the surface and for the production rates of molecules. In the first stage of 
the experiment the fluxes were Fh = -Fb, while the diffusion and desorption rates were 
practically zero, due to the low surface temperature. In the second stage the two beams 
were turned off, while the diffusion and desorption coefficients quickly increased due to the 
heating of the surface. 

The isotope effect studied here gives rise to a small difference in the energy barriers 
for diffusion of H and D atoms, -E'o(H) < Eq{D), and a larger difference in the desorption 
energies, where -Ei(H) < Ei{D). We will now examine how these differences would affect the 



32 A. Lipshtat, 0. Biham and E. Herbst 

results of the TPD experiment. We first consider the effect of the difference in the diffusion 
barriers. The rate of formation of HD molecules is described by Eq. 0, where the diffusion 
rates of the H and D atoms sum up. It is thus clear that even if the D atoms are immobile, 
this gives rise only to about a factor of 2 in the production rate compared to the case in which 
Ad = An. Such differences have little effect on the energy barrier for diffusion, Eq, obtained 
from fitting the experimental TPD curve. This is because Eq appears in the exponent of 
Eq. Therefore, when the experimental TPD curves are fitted by a rate equation model 
that includes a single energy barrier for diffusion, Eq, the result is very close to the lower 
diffusion barrier among the H and D, namely to -E'o(H). 

Consider the case in which the desorption barrier of D atoms, Ei{D), is larger than the 
desorption barrier of H atoms, i?i(H). In this C3iSG, clS the temperature is raised during a 
TPD experiment, D atoms desorb more slowly than H atoms. The formation rate of HD 
molecules is thus limited by the availability of H atoms. When the resulting TPD curves are 
fitted by a rate equation model that includes only a single energy barrier for desorption, Ei, 
the fitted value turns out to be a very good approximation to -E'i(II). 

It thus turns out that slower diffusion and desorption rates of D atoms vs. H atoms, have 
little effect on the TPD curves that describe the rate of formation of HD molecules. When 
these TPD curves are fitted by a rate equation model with a single diffusion barrier, Eq, 
and a single desorption barrier, Ei, these barriers provide good approximations to -Eo(H) 
and £'i(H). This conclusion is not limited to the case in which the experiment is done on a 
macroscopic surface. It is expected to be valid even if the experiment is done on a powder 
that consists of small grains. 

The origin of the isotope effect on dust grains in the interstellar medium is the long 
residence time of D atoms on the grain compared with H atoms. As a result, an (H or D) 
atom that hits a grain is more likely to encounter a D atoms on the grain than expected 
based on the gas phase abundance ratio. In the TPD experiments all the H and D atoms are 
deposited on the surface at the initial stage of the experiment, when the surface temperature 
is low and the adsorbed atoms are immobile. During the heating stage, when the adsorbed 
atoms become mobile, the beams are turned off and there is no fresh supply of H or D atoms. 
Therefore, the D atoms that remain on the surface after most of the H atoms are desorbed 
tend to form D2 molecules. The isotope effect is thus expected to affect the TPD curve of 
D2 molecules. However, it does not affect the TPD curve of HD molecules, which is the one 
that was measured in the experiments. 



